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Abstract 

The basic intent of this paper is to model 70 Ophiuchi A using the latest asteroseismic 
observations as complementary constraints and to determine the fundamental parameters 
of the star. Additionally, we propose a new quantity to lift the degeneracy between the initial 
chemical composition and stellar age. Using the Yale stellar evolution code (YREC7), we 
construct a series of stellar evolutionary tracks for the mass range M = 0.85 - 0.93 M Q with 
different composition (0.26 - 0.30) and Z, (0.017 - 0.023). Along these tracks, we select a 
grid of stellar model candidates that fall within the error box in the HR diagram to calculate 
the theoretical frequencies, the large- and small- frequency separations using the Guenther's 
stellar pulsation code. Following the asymptotic formula of stellar /?-modes, we define a 
quantity roi which is correlated with stellar age. Also, we test it by theoretical adiabatic 
frequencies of many models. Many detailed models of 70 Ophiuchi A have been listed 
in Table 3. By combining all non-asteroseismic observations available for 70 Ophiuchi A 
with these seismological data, we think that Model 60, Model 125 and Model 126, listed 
in Table 3, are the optimum models presently. Meanwhile, we predict that the radius of this 
star is about 0.860 - 0.865 R Q and the age is about 6.8 - 7.0 Gyr with mass 0.89 - 0.90 M Q . 
Additionally, we prove that the new quantity roi can be a useful indicator of stellar age. 
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1 Introduction 

The solar five-minute oscillations have led to a wealth of information about the in- 
ternal structure of the Sun. These results stimulated various attempts to detect solar- 
like oscillations for a handful of solar-type stars. Individual p-mode frequencies 
have been identified for a few stars: a Cen A (Bouchy and Carrier, 2002; Bedding 
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et al., 2004), a Cen B (Carrier and Bourban, 2003a; Kjeldsen et al., 2005), /i Arae 
(Bouchy et al., 2005), HD 49933 (Mosser et al., 2005), fi Vir (Martic et al, 2004a; 
Carrier et al., 2005b), Procyon A (Martic et al, 2004b; Eggenberger et al., 2004a), 
77 Bootis (Kjeldsen et al., 2003; Carrier et al., 2005a), /? Hyi (Bedding et al, 2001; 
Carrier et al., 2001) and 6 Eri (Carrier et al., 2003b). Based on these asteroseismic 
data, numerous theoretical analyses have been performed in order to determine pre- 
cise global stellar parameters and to test the various complicating physical effects 
on the stellar structure and evolutionary theory (Thevenin et al., 2002; Eggenberger 
et al., 2004b, 2005; Kervella et al., 2004; Miglio and Montalban, 2005; Provost et 
al., 2004, 2006). 

Recently, Carrier and Eggenberger (2006) detected solar-like oscillations on the 
K0 V star 70 Ophiuchi A (HD 165341), and identified some possible existing fre- 
quencies. They obtained the large separation Av = 161.7 + 0.3 jxHz by observation 
over 6 nights with HARPS. The spectroscopic visual binary system 70 Ophiuchi 
is one of our nearest neighbors (5 pc) and is among the first discovered binary 
stars. It was observed first by Herschel in 1779. So 70 Ophiuchi A is famous as 
the primary of a visual and spectroscopic binary system in the solar neighborhood. 
Although many observation data have been obtained since 1779, the theoretical 
analysis of 70 Ophiuchi A has only been made by Fernandes et al. (1998). By a 
calibration method which take into account the simultaneous evolution of the two 
members of the binary system, they analyzed the 70 Ophiuchi A by means of stan- 
dard evolutionary stellar models using the CESAM code bf (Morel, 1997) without 
microscopic diffusion. They found that the metallicity of 70 Ophiuchi A is very 
close to the solar one, the values of mixing-length parameter a and helium abun- 
dance Y are near the Sun. They thought that the star is younger than the Sun and 
3 + 2 Gyr is probably an limit considering the age versus stellar rotation relation 
with its rotation velocity (vsini « 16 km s~ l ). 

The aim of our paper is to present the model which can be constrained by these 
seismology data. The observational constraints available for 70 Ophiuchi A are 
summarized in Sect. 2, while the numerical calculations are presented in Sect. 3. 
The seismic analyses are carried out and a new quantity r i as a indication of stellar 
age is proposed in Sect. 4. Finally, the discussion and conclusions are given in Sect. 
5. 



2 Observational Constraints 

2.1 Non-asteroseismic observation constraints 

The mass of this star was investigated by Batten et al. (1984), Heintz et al. (1988), 
Fernandes et al. (1998) and Pourbaix et al. (2000), respectively. In the paper, we 
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Table 1 



Non-asteroseismic observational data of 70 Ophiuchi A. 



Observable 


Value 


Source 


Mass M/M Q 


0.89 ± 0.04 


(1) 


Effective temperature T e ff(K) 


5322 ± 20 


(2) 


Luminosity log(L/L Q ) 


-0.29 ± 0.03 


(1) 


Metallicity [Fe/H] s 


0.0 ±0.1 


(1) 


Surface heavy element abundance [Z/X] s 


0.02365 ± 0.00535 


(3) 



References. — (1) Fernandes et al. (1998), (2) Gray and Johnson (1991), (3) this paper. 

adopt the value of mass deduced from Fernandes et al. (1998). The effective tem- 
perature was determined by Gray and Johnson (1991). So far, the metallicity ob- 
tained by observation are [Fe/H] = -0.05 (Peterson, 1978) and [Fe/H] = 0.00 (Perrin 
et al., 1975). We choose the [Fe/H] = 0.0 + 0.1 as a representative value according 
to Fernandes et al. (1998). 

The mass fraction of heavy elements, Z, was derived assuming log[Z/X] « [Fe/H] + 
log[Z/X] Q and [Z/X] Q = 0.0230 (Grevesse and Sauval, 1998), for the solar mixture. 
So we can deduce the [Z/X] s = 0.0183 - 0.0290. 

All non-asteroseismic observational constraints are listed in Table 1 . 



2.2 Asteroseismic constraints 



Solar-like oscillations of 70 Ophiuchi A have been detected by Carrier and Eggen- 
berger (2006) with the HARPS spectrograph. Fourteen individual modes are identi- 
fied with amplitudes in the range 1 1 to 14 cm s~ l . Although they listed two groups 
of frequencies by mode identification (see Table 2 in Carrier and Eggenberger, 
2006), one group of frequencies with a average large separation Av = 161.7 jxHz 
was suggested to be more reliable than the other with a average large separation 
Ay = 172.2 ixHz. The star 70 Ophiuchi A is very similar to a Cen B with the same 
spectral type and similar large separation, which has a mean small separation of 
10 juHz. It is thought that the small separation should be similar. By inspecting 
the results of the mode identification, they note that the value of the small sep- 
aration coming from the identification with the large separation of 172.2 /j.Hz is 
significantly different from 10 [iHz. If the large separation is 172.2 /uHz, the small 
separation will be lower than 6.5 /uHz in the frequency range 3 - 4.5 mHz. Although 
this identification is less reliable than the one with a large separation of 161.7 f^Hz, 
the solution Av = 172.2 fiHz can not be ruled out definitely. We refer to these two 
groups of results in the paper and make analyses in Sect. 4 and Sect. 5. 
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Table 2 



Input parameters for model tracks. 



Variable 


Minimum Value 


Maximum Value 


6 


Mass M/M Q 


0.85 


0.93 


0.01 


Initial heavy element abundance Z, 


0.017 


0.023 


0.001 


Initial Helium abundance F, 


0.26 


0.30 


0.01 



Note. — The value 6 defines the increment between minimum and maximum parameter val- 
ues used to create the model array. 




3.735 3730 3.725 3.720 2 4 6 8 10 12 

Log(Teff) Age [Gyr] 



Fig. 1. Evolutionary tracks in the HR diagram and^ 2 as a function of age for 70 Ophiuchi 
A. a) The selected individual stellar evolutionary tracks (44 in total). b)x 2 values calculated 
for 70 Ophiuchi A observational data using different Z„ F, and Mass, plotted as a function 
of age. x 2 refers to non-asteroseismological observables as denoted by eq. (1). 

3 Stellar models 

We will construct a grid of stellar evolutionary models by Yale stellar evolution 
code (YREC; Guenther et al., 1992) with microscopic diffusion. The initial zero- 
age main sequence (ZAMS) model used for 70 Ophiuchi A was created from pre- 
main sequence evolution calculations. In these computations, we do not consider 
rotation and magnetic field effect. These models are computed using OPAL equa- 
tion of state tables EOS 2001 (Rogers and Nayfonov, 2002), the opacities interpo- 
lated between OPAL GN93 (Iglesias and Rogers, 1996) and low temperature tables 
(Alexander and Ferguson, 1994). Using the standard mixing-length theory, we set 
a = 1.7 for all models, close to the value which is required to reproduce the solar 
radius under the same physical assumptions and stellar evolution code. Meanwhile, 
it must be emphasized that there are still a number of uncertainties in our analyses, 
foremost among which is the still open question of mixing-length theory responsi- 
ble for the stellar model. The nuclear reaction rates have been updated according 



4 



— , 1 1 1 1 1 1 1 r 

(b) Model 60 




3000 4500 

V [fiUz] 



6000 3000 



4500 



6000 



Fig. 2. Large-frequency separations vs. frequency for the Models 54, 60, 125, 126 (in Table 
3). The observable large separations Ay versus frequency for /?-modes of degree I = 0(»), 
/ = l(A)and I = 2(b) are obtained from Carrier and Eggenberger (2006), which correspond 
to the average large separation of 161.7 fiHz- Open symbols correspond to large separation 
averages taken between non-successive modes and vice versa. All individual errors are 
fixed to V2 x 1.1 fiHz (half resolution). 



to Bahcall and Pinsonneault (1995). The Krishna-Swamy Atmosphere T-r relation 
is used for this solar- like star (Guenther and Demarque, 2000). Also, we consider 
the microscopic diffusion effect, by using the diffusion coefficients of Thoul et al. 
(1994). Since 70 Ophiuchi A, like a Cen B, is less massive than the Sun, the mass 
contained in its convective zone is much larger and, therefore, the effect of micro- 
scopic diffusion is much smaller (Miglio and Montalban, 2005; Morel and Baglin, 
1999). However, it is necessary to consider this effect as physical process in stellar 
modeling (see Provost et al., 2005, 2006). 



In general, the determination of parameters (M, t, Y t , Z,) fitting the observational 
constraints needs two steps. The first step is to construct a grid of models with 
position in the HR diagram in agreement with the observational values of the lu- 
minosity, the effective temperature and the surface metallicity. The principal con- 
straints deduced from non-asteroseismic observation are listed in Table 1 . The error 
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box, which is composed of observational effective temperature and luminosity, rep- 
resents the possible position of 70 Ophiuchi A in the HR diagram (see Fig. la). 
According to the results of Fernandes et al. (1998), we list the parameter space of 
mass M, the initial heavy-element abundance Z, and the initial helium abundance 
Yi in Table 2. Since the microscopic diffusion is included in our paper, we give the 
wider parameter space of initial heavy-element abundance Z, than the range of Z, 
of Fernandes et al. (1998). 

By adjusting three parameters M, Y t and Z, listed in Table 2, we can obtain many 
evolutionary tracks passing through the error box in the HR diagram. Now we con- 
sider a function which describes the agreement between the observations and the 
theoretical results: 



where C represent the following quantities: L/L Q , T e ff and [Z/X] s , C represent 
the theoretical values and C obs represent the observational values listed in Table 1 . 
The vector crC" bs contains the errors on these observations which are also given in 



As Fernandes(1998) has pointed that the age of 70 Ophiuchi A is 3 + 2 Gyr, it is 
reasonable for us to choose the evolutionary tracks passing through the error box 
within 8 Gyr. We select 44 evolutionary tracks passing through error box as our 
possible candidates to go on with our investigations. Fig. la gives 44 evolutionary 
tracks, and Fig. lb presents^- 2 as a function of the age correspondingly. It is well- 
known that^- 2 is smaller, the more competitive is the candidate. Fig. lb shows that 
models with^ 2 smaller than 1 have ages between 3 Gyr to 7Gyr. From Fig. la, we 
find that the upper section of error box is empty. The reason of the empty upper 
section of the error box is related to the range of initial parameters, like mass, 
initial composition and specially the mixing length parameter. We think that the 
future interferometric measurement of the radius could reduce the domain of the 
possible position in the HR diagram (e.g., Provost et al., 2006). 

The second step is to determine the optimum model using the asteroseismic mea- 
surements. We will select a grid of models along these 44 tracks shown in Fig. 
la to calculate the low-/ /?-modes frequencies. We list the representative models 
extracted from every tracks in Table 3. 

The detailed pulsation analysis is described in the next section. 
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Fig. 3. Differences between calculated and observable frequencies for the Model 54, 60, 
125, 126 in Table 3. The systematic shifts (D v ) for the four models are 30.059, 63.75, 
70.252, 63.24 [iHz„ respectively . The observable frequencies correspond to the average 
large separation of 161.7 \iHz (see text for details. 

4 Pulsation analysis 



4. 1 Selecting the optimum model 



Using Guenther's pulsation code (Guenther, 1994), we calculate the adiabatic low-/ 
jp-mode frequencies of the selected models. We define the large separations Av and 
small separations Sv in the usual way (Tassoul, 1980): 

Av n>/ = v nJ - v n _ u (2) 

and 

5v nJ = V n<t - V n -ij + 2, (3) 

where n is the radial order, / is the degree, and v is the frequency. Because the 
expected acoustic cutoff has a limit, we only calculate the mean large- and small- 
separations by averaging over n = 10-30 (See Murphy et al., 2004). Within these 



7 



6000 
5500 

- 5000 

N 

5. 4500 

* 4000 
3500 

§888 

5500 
^ 5000 

N 

t 4500 

* 4000 

3500 
3000 



-s— i- 



-i— i- 



(c) 



A' A 
A A 

A 
A 
A 
A A 

A 

A 

A 

A 
ZA 
A 
A 
A 
A 
A 

A 

■ A| , 



A 
A 

A 

A 

A 
A 
A 
A 
A 
A 
A 

A 

A 

A 
A 
A 
A 

H — H 1 - 



Model 



-i— i— |— r 

lodel 54 



□ o 

□ 

□ 0- 

□ 

□ 

□ ■ 

□ 

□ 

□ ■ 

□ 

□ 

□ ■ 

□ 

□ 

□ 0— 

□ o: 
, □ - 



I I t 

1(20 DO 

mo 

ao 

□€> 

□ t 

□ o 
mo 

□ o 

□ OJ 

□ ot 

□ o 

□ o • 

□ o 

□ • 

□ Ot 

b a i 



-I— I- 



(b) 



7T 

A A 

A 
A 
A 
AA 

A 

A 

A 
A 
ZA 
A 
A 
A 
A 

A 

A 



T 



Model 60 



-l— l- 



-l— l- 



(d) 



A A 
A A 

A 
A 
A 
AA 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 



H — h 



□ b 

□ o 

□ 
■ 

□e 
□c 

□ • 

□ o 

MO 

□ O 

□ » 

□ 0* 

□ O 

□ 09 

□ O 

□ ot 

□ cm 



-\—\- 



□ 

Model 126 no 
■ o 

□e 

□ • 

□ o 
no 

□ o 

□ at 

□ ot 

□ o 

□ ot 

□ 

□ Ot 

□ ot 

I , , I p p I 



• l(obs)=0 
i l(obs)=1 
■ l(obs)=2 
o l(theo)-0 
a l(theo) = 1 
□ l(theo)=2 



30 60 90 120 150 
v modulo kv [//Hz] 



30 60 90 120 150 
v modulo kv [//Hz] 



Fig. 4. Echelle diagrams for the Models 54, 60, 125, 126 (in Table 3), which average large 
separations (Av) are 166.22, 161.7, 161.92, 161.68 fiHz, respectively. The systematic shift 
(D v ), which have been applied to the theoretical frequencies, are 30.059, 63.75, 70.252, 
63.24 /uHz, for the four models respectively. Open symbols refer to the theoretical frequen- 
cies, and filled symbols to the observable frequencies. Circles are used for I - modes, 
triangles for / = 1 modes, and squares for / = 2 modes. The observable frequencies corre- 
spond to the average large separation of 161.7 /uHz (see text for details). 



44 tracks, we list 129 models in Table 3. (Av/) represents the mean of large sep- 
arations Av n j for n = 10 to 30. The frequency range corresponds to about 2000 
fiHZ-6000 /uHz. Additionally, (Av) represents the mean of (Av/) for / = to 3. 
In the same way, (Svo2> and (t>vi3> represent the mean of £v„,o and 6v n j for n = 
10 to 30, respectively. So far, we only know the large separations and the four- 
teen individual modes of the star based on the asteroseismic data of Carrier and 
Eggenberger (2006). Guenther (1998) pointed that the large separations are most 
easily identifiable characteristics in the p-mode spectrum. Because they are seen 
as a peak in the Fourier transform of the power spectrum and they are mostly un- 
contaminated by composition effects, these large separations provide an efficient 
way to constrain stellar model. It is also important to remember that the theoretical 
frequencies calculated in our paper should not be expected to match the observed 
frequencies of Carrier and Eggenberger (2006) perfectly. We think that there are 
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three reasons. Firstly, our theoretical models do not match the mass and radius of 
the star precisely. Secondly, the uncertainty in calculating the sound speed in the 
outer layers of the models comes into being, where non-adiabatic effects become 
important. Thirdly, at high frequencies, the effect of the convection-oscillation in- 
teractions is larger and the description of convection is open problem. Although 
the differences between the theoretical frequencies and the observed frequencies 
could result in significant effect on the large separations, we think that the effect is 
small due to the large separations correspond to differences between frequencies of 
modes with the same angular degree / and consecutive radial order n. Therefore, in 
our paper, we think that the matching the observable large separations is the impor- 
tant criterion to select the optimum model. In Table 3, we find that the average large 
separations of Model 60, Model 125 and Model 126 are 161.7, 161.92 and 161.68 
liHz, in good agreement with the mean value derived from Carrier and Eggenberger 
(2006). So we can tentatively say that these models may be the best fit models. In 
Fig. 2, we plot the observational results about the large separations and the errors. 
Also we plot the large separation as a function of frequency for Model 54 in Fig. 
2a, Model 60 in Fig. 2b, Model 125 in Fig. 2c and Model 126 in Fig. 2d. We clearly 
find that the theoretical large separations of the Model 60, Model 125 and Model 
126 are consistent with the observations. Model 54, as the representative of many 
non-fit models, is not consistent with the observational large separations. There- 
fore, we have sufficient reasons to say that Model 60 , Model 125 and Model 126 
are really the best fit models. Meanwhile, we can predict that the radius of star is 
0.860 - 0.865 R Q and the age is about 6.8 - 7.0 Gyr with mass 0.89 - 0.90 M 
presently. 

Once the asteroseismic observation can confirm the large separations to be 161.7 + 
0.3 ixHz and the theory Model 60 , Model 125, Model 126 are considered as the best 
models, we can predict that the mean small separation (6vo 2 ) is about 10.29 - 10.48 
jxHz and the radius of the star is about 0.860 - 0.865 R Q . Direct measurements of 
stellar diameters from interferometric observations should provide an independent 
check for asteroseismic predictions such as Kervella et al. (2003a, 2003b). 

In order to compare the theoretical p-mode frequencies deduced from the models 
in Table 3 with the observational frequencies provided by Carrier and Eggenber 
(2006), we plot the echelle diagram of every model and find that no model can 
fit observational frequencies. For the exact values of the frequencies, considering 
above three reasons, a linear shift of a few fiHz between theoretical and obser- 
vational frequencies is perfectly acceptable. Taking into account it, we define the 
mean value of the difference between the theoretical and observational frequencies 
(e.g., Eggenberger et al., 2004b, 2005): 




(4) 



where N is the number of observable frequencies (N = 14). 
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Fig. 5. The "reduced" large separation vs. mass for each of the 129 stellar models. 

Takeing into account the systematic difference (D v ) between theoretical and ob- 
servable frequencies, we plot the differences between calculated and observed fre- 
quencies in Fig. 3 and the echelle diagram in Fig. 4. The observable frequencies 
correspond to the average large separation of 161.7 [iHz in these figures. Fig. 3a, 
Fig. 3b, Fig. 3c and Fig. 3d correspond to the Model 54 with (D v ) = 30.059 /uHz, 
Model 60 with (D v ) = 63.75 fiHz, Model 125 with (D v ) = 70.252 /uHz, and Model 
126 with (D v ) = 63.24 pHz, respectively. Fig. 4a, Fig. 4b, Fig. 4c and Fig. 4d show 
the echelle diagram of the Model 54, Model 60, Model 125 and Model 126 re- 
spectively. For p-modes in the asymptotic theory (n » /), the large separations are 
nearly constant; meanwhile the so-called "echelle diagrams" present the frequen- 
cies in ordinates, and the same frequencies modulo the average large separation in 
abscissa. So the asymptotic theory predicts an approximated vertical line for given 
degree. In this case, Fig. 4b, Fig. 4c and Fig. 4d show that the theoretical frequen- 
cies of Model 60, Model 125, Model 126 can fit the observable frequencies with 
161.7 ixHz. Meanwhile, we find that the systematic differences (D v ) are larger than 
the results of a Cen B obtained by Eggenberger (2004b). It is interesting to analyze 
the difference in future work. 



4.2 Asymptotic formula and frequency analysis 



4.2.1 Large Separations and Small Separations 

It is well known from asymptotic theory that the large separations are mainly sen- 
sitive to the stellar radius (Tassoul, 1980; Christensen-Dalsgaard, 1984). More pre- 
cisely, the asymptotic behavior of Av is expected to scale with (M/R 3 ) l/2 , where M 
is the mass of the star and R is its radius. Meanwhile, Murphy et al. (2004) find 
that a degeneracy in predicted radius occurs for models of different mass. Here, the 
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Fig. 6. The ratio of small separations adjacent in / vs. age for each of 129 stellar models. 

degeneracy means that the (Av) changes with radius and mass (see Fig. 4 in Mur- 
phy et al., 2004). In order to lift the degeneracy, Fernandes and Monteiro (2003) 
and Murphy et al. (2004) assumed homology to compare theoretical models by 
introducing a "reduced" radius, such as 

(Av r ) = (Av n MR/Ref 2 . (5) 

Here, we name the quantity (Av r ) "reduced" large separation. We draw the "re- 
duced" large separation (Av r ) versus mass in Fig. 5 and list the values of (Av r ) for 
each model in Table 3. From Fig. 5, we find that the degeneracy was lifted approx- 
imately. It is easily seen that the values of the (Av r ) are relatively consistent with 
each mass. It is successful using the (Av r ) instead of large separations to indicate 
the stellar mass. 

The small separations, like the large separations , will be visible as peaks in the 
Fourier transform of the power spectrum. At the earlier stage, Christensen-Dalsgaard 
(1984) proposed that the calculation of small separations could put a constraint on 
the age the star. Subsequently, Ulrich (1986) realized that only if the composition of 
the star is known completely can one use the small separations to correctly identify 
a stellar age. This point has been illustrated in Murphy et al. (2004). Thus, the var- 
ious chemical compositions create a degeneracy in age determination (see Murphy 
et al., 2004). Namely, the small separations Sv change with the initial composition 
and age. In the next section, we will discuss this problem and propose a quantity 
which may be correlated with stellar age. 

4.2.2 A new quantity be proposed 

At the present time, we can know the stellar internal structure and understand the 
stellar evolution from oscillation frequencies. Thus asteroseismology provides a 
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window to "see" the interior of star. But the observation of solar-like oscillation is 
very difficult because of their small amplitude. So far, we only obtain the knowledge 
of the stellar interior from the limited modes (/ = 0, 1, 2, 3) which can be observed. 
Many authors proposed some quantity as diagnostic purposes to probe the stellar 
internal and constraint the model parameters (Christensen-Dalsgaard, 1984, 1988, 
1993; Ulrich, 1986, 1988; Gough, 1987, 1990a, 2003). For p-modes of solar-like 
stars, the usual frequency separations are the large separation defined by equation 
(2) and the small separation defined by equation (3). Additionally, Roxburgh (1993) 
and Roxburgh and Vorontsov (2003) considered the following separations: 

8oi(n) = \{v„-\ fi - 4v ra _ u + 6v„, - 4v n>1 + v, 1+U) ), (6) 



Sio(n) = -^(Vn-i.i - 4v ;i _i j0 + 6v nJ - 4v n+1>0 + v n+l>l ), (7) 
and defined the ratios dy of small to large separations as follows: 



doiin) = - — — , d l3 (n) = 



Ayi(n)' 1JV Av (n + 1)' 

6v m (n) 6v l0 (n) 
"oi in) = - — — , d w {n) = 



(8) 



Av^nY Av (/z + 1)" 



The ratios d {j of small to large separations are independent of the structure of the 
outer layers of a star, and therefore provide a diagnostic of the stellar interior alone. 

In addition, Gough (1990a), Monteiro and Thompson (1998), Vauclair and Theado 
(2004), Houdek and Gough (2007a) gave the second differences A 2 v ; (n): 



A 2 vi(n) = v„ +u + v„_ u - 2v nJ . (9) 

The second differences A 2 v/(n) can be used to reveal the variation of the first adi- 
abatic exponent yi dependent of the influence of the ionization of helium on the 
low-degree acoustic oscillation frequencies in model of solar-type stars. Recently, 
Houdek and Gough (2007b) stated that the second differences can provide a mea- 
sure of helium abundance and hence precisely lift the degeneracy between compo- 
sition and age. 

Summarizing the above character separation, we find that investigation of lifting 
the degeneracy between the chemical compositions and the age is interesting. We 
begin with our investigation from a well-known asymptotic formula. 

The asymptotic formula for the frequency v n j of a stellar p-mode of order n and 
degree / was given by Tassoul (1980): 

v n ,i - in + l - + 6)v () - [Al(l + 1) - fflvfoj, (10) 
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where, the characteristic vo is related to the run of sound travel time across the 
stellar diameter; A is a measure of the sound-speed gradient and most sensitive 
to conditions in the stellar core (see Gough and Novotny, 1990b; Gough, 2003; 
Christensen-Dalsgaard, 1993; Guenther and Brown, 2004), e and B are constants 
which are the functions of the equilibrium model. It should be noted that the classi- 
cal asymptotic theory of Tassoul (1980), although providing good results at the first 
order in frequency, does not represent with accuracy the p-mode spectrum of the 
stars considered. Several authors (e.g., Gabriel, 1989; Audard and Provost, 1994; 
Roxburgh and Vorontsov, 2000a, 2000b, 2001) have been discussed the difficulties 
of the asymptotic theory, particularly for evolved models with rapid variation in the 
sound speed in the core. 

Using the equation (2) and the asymptotic formula (10), the large separation can be 
written as follows (Gough and Novotny, 1990b): 

Av„,, = V n j - V n -u 

r ai,i i\ m 2 V n-U~ v n,l (11) 

= Vo - [Al(l + 1) - #K , 

y n ,i ■ v «-u 

Taking the first order of v n j for the n » / approximately, we can obtain the result 
like Gough and Novotny (1990b) and equation (11) becomes:, 



Aig+p-B ! 

Av„,/ w v [l - ■ -] 

(n + i + e)(n - 1 + { + e) (12) 

« V(). 

Using the same approximate method, we can obtain the expression of small sepa- 
ration: 



5v n / = 6v, i +2 



(41 + 6)v A 



n + j + e 



(13) 



Due to the small separations are rather sensitive to composition and therefore to 
the structure of the core, especially the extreme sensitivity of the stellar core den- 
sity stratification to several parameters (Guenther and Demarque, 2000; Morel et 
al., 2000), we define another quantity about the ratio of average small separation 
adjacent in /: 



ro, = (14) 

Using the equation (14), we calculate the values of r m and list it in table 3. Based 
on the results of numerical calculations, we plot the ratio r i versus age in Fig. 6. 
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Fortunately, we find that the ratio r m is tightly correlated with age and decreases 
monotonously with age. We think that the likely reason comes from the perturba- 
tion to the gravitational potential, neglected in the asymptotic relation (10), which 
affects modes of the lowest degrees most strongly and which probably increases 
with evolution due to the increasing central density. These effects are most impor- 
tant for modes of the lowest degrees which penetrate most deeply and hence affect 
£v ,2 more than 6v\$, leading to the dependence of r i on age. 

From the Fig. 6, the values of r 0; i in table 3 and the above discussion, we can 
conclude that this quantity r 0l can lift the degeneracy between the chemical com- 
positions and age. The analysis was inspired by Fernandes and Monteiro (2003) 
and Murphy et al. (2004). So, we can obtain the r i which may indicate stellar 
age, if we consider a frequency ratio. As illustrated in Fig. 6, the quantity r m is 
tightly correlate with stellar age over a substantial range of the remaining parame- 
ters, including composition. At the same time, we need to point out that the range 
of variation of this quantity is relatively modest, compared to likely observational 
errors. Also, it is clear that the present data for 70 Ophiuchi A are not adequate to 
evaluate this quantity. We think that using this quantity to evaluate the stellar age 
will be convenient based on asteroseismic data which will be provided in the future. 



5 Discussion and conclusions 

The models of 70 Ophiuchi A are obtained by fitting effective temperature, lumi- 
nosity, surface metallicity and asteroseismic observations. 

We list a series of possible models in Table 3. So far, we think that Model 60, 
Model 125 and Model 126 are the best models. With the advance in observation, 
the precise asteroseismic data will provide more strict constraints on theoretical 
models. The conclusions of the paper are: 

1 . Using the latest asteroseismic observation, we try our best to construct the best 
model of 70 Ophiuchi A. So far, we only select the Model 60, Model 125 and Model 
126, which can be fit for the observation, as the optimum models. These models are 
correspond to the radius of this star is about 0.860 - 0.865 R Q and the age is about 
6.8 - 7.0 Gyr with mass 0.89 - 0.90 M . 

2. By calculating many theoretical models, we want to use the theoretical frequen- 
cies to compare with observational frequencies and help to definitely validate the 
observational data. 

3. We obtain a new quantity r m which can lift the degeneracy between the initial 
compositions and the stellar age. By calculation, we prove that it can be valuable 
for the indication of the stellar age. 
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4. Important point is that we test our stellar structure and evolution theory. Mean- 
while we find the confrontation of observations and theoretical models. Actually 
we have neglect some important effects, such as rotation and magnetic field, which 
can impact on the internal structure and evolution. Thus, the theory models which 
we have constructed can not fit observation perfectly. The detailed comparisons of 
individual mode frequencies will also require taking into account the effects of tur- 
bulence in the outer convective unstable layers in the stellar models, which shift 
the observed frequencies. The parameterization of turbulence tested on the Sun by 
Li et al. (2002) can be applied to models for solar-like stars as well. This param- 
eterization can be extended to other stars by using the three-dimensional radiative 
hydrodynamic simulations of Robinson et al. (2003), which are based on the same 
microscopic physics and can readily be parameterized in the YREC stellar evolu- 
tion code. 
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Table 3 

Model parameters (M: mass; Yf. initial Helium abundance; Z ; : initial heavy element 
abundance; \ogT e ff log-effective temperature; logL: log-luminosity; R: radius; (Ay/>: 
Ti 2n°=io Av nf, (Av): iS? =0 <Av / >; <<5v 02 >: £ Zf =10 Sv nfl ; <<5v 13 ): £ Z* 10 <5v„,i ; r 01 : fg). 
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Table 3 
-Continued 
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128.9844 


0.6020 


62 


0.89 


0.27 


0.021 


3.7276 


-0.309 


0.8195 


3.45 


173.71 


173.74 


173.89 


174.21 


173.88 


13.65 


22.80 


128.9951 


0.5987 


63 


0.90 


0.26 


0.019 


3.7273 


-0.319 


0.8118 


1.90 


177.12 


177.13 


177.29 


177.62 


177.29 


15.07 


24.75 


129.6755 


0.6089 


64 


0.90 


0.26 


0.019 


3.7276 


-0.315 


0.8144 


2.10 


176.24 


176.26 


176.42 


176.74 


176.41 


14.87 


24.47 


129.6522 


0.6077 


65 


0.90 


0.26 


0.020 


3.7261 


-0.319 


0.8159 


2.55 


175.82 


175.84 


175.99 


176.32 


175.99 


14.47 


23.92 


129.7010 


0.6049 


66 


0.90 


0.26 


0.020 


3.7270 


-0.308 


0.8233 


3.10 


173.44 


173.47 


173.62 


173.94 


173.61 


13.92 


23.15 


129.6916 


0.6013 


67 


0.90 


0.26 


0.020 


3.7276 


-0.300 


0.8282 


3.45 


171.90 


171.93 


172.08 


172.40 


172.07 


13.56 


22.66 


129.6904 


0.5984 


68 


0.90 


0.26 


0.021 


3.7250 


-0.319 


0.8203 


3.25 


174.43 


174.46 


174.61 


174.93 


174.60 


13.85 


23.05 


129.7189 


0.6009 


69 


0.90 


0.26 


0.021 


3.7254 


-0.314 


0.8237 


3.50 


173.37 


173.39 


173.54 


173.87 


173.54 


13.60 


22.71 


129.7338 


0.5989 


70 


0.90 


0.26 


0.021 


3.7259 


-0.308 


0.8271 


3.75 


172.29 


172.32 


172.46 


172.79 


172.46 


13.35 


22.36 


129.7255 


0.5970 


71 


0.90 


0.26 


0.021 


3.7267 


-0.298 


0.8342 


4.25 


170.09 


170.13 


170.27 


170.60 


170.27 


12.84 


21.66 


129.7309 


0.5928 


72 


0.90 


0.26 


0.021 


3.7276 


-0.285 


0.8424 


4.80 


167.62 


167.67 


167.81 


168.13 


167.80 


12.28 


20.87 


129.7387 


0.5884 


73 


0.91 


0.26 


0.020 


3.7274 


-0.319 


0.8112 


1.40 


178.34 


178.35 


178.51 


178.84 


178.51 


15.46 


25.28 


130.4231 


0.6116 


74 


0.91 


0.26 


0.020 


3.7277 


-0.315 


0.8138 


1.60 


177.46 


177.47 


177.63 


177.95 


177.62 


15.26 


25.00 


130.3972 


0.6104 


75 


0.91 


0.26 


0.021 


3.7262 


-0.319 


0.8154 


2.05 


176.98 


176.99 


177.15 


177.47 


177.14 


14.85 


24.43 


130.4286 


0.6079 


76 


0.91 


0.26 


0.021 


3.7270 


-0.309 


0.8221 


2.55 


174.82 


174.84 


174.99 


175.31 


174.99 


14.36 


23.74 


130.4368 


0.6049 


77 


0.91 


0.26 


0.021 


3.7277 


-0.300 


0.8277 


2.95 


173.07 


173.09 


173.24 


173.56 


173.24 


13.95 


23.17 


130.4541 


0.6021 


78 


0.85 


0.29 


0.022 


3.7244 


-0.353 


0.7914 


4.05 


178.88 


178.91 


179.06 


179.39 


179.06 


13.56 


22.78 


126.0645 


0.5953 


79 


0.85 


0.29 


0.022 


3.7261 


-0.332 


0.8045 


5.00 


174.56 


174.62 


174.75 


175.08 


174.75 


12.59 


21.42 


126.0974 


0.5878 


80 


0.85 


0.29 


0.022 


3.7277 


-0.310 


0.8195 


6.00 


169.78 


169.87 


169.99 


170.32 


169.99 


11.52 


19.94 


126.1092 


0.5777 


81 


0.86 


0.28 


0.022 


3.7248 


-0.323 


0.8176 


5.70 


171.40 


171.47 


171.60 


171.93 


171.60 


11.91 


20.46 


126.8611 


0.5821 


82 


0.86 


0.28 


0.022 


3.7262 


-0.305 


0.8297 


6.50 


167.68 


167.77 


167.89 


168.22 


167.89 


11.08 


19.31 


126.8839 


0.5738 


83 


0.86 


0.28 


0.022 


3.7277 


-0.284 


0.8446 


7.40 


163.28 


163.40 


163.51 


163.85 


163.51 


10.13 


17.97 


126.9173 


0.5637 


84 


0.87 


0.28 


0.022 


3.7263 


-0.321 


0.8141 


4.55 


173.48 


173.53 


173.67 


174.00 


173.67 


12.82 


21.70 


127.5679 


0.5908 


85 


0.87 


0.28 


0.022 


3.7270 


-0.312 


0.8199 


4.95 


171.66 


171.72 


171.85 


172.18 


171.85 


12.41 


21.13 


127.5824 


0.5873 


86 


0.87 


0.28 


0.022 


3.7277 


-0.303 


0.8260 


5.35 


169.78 


169.84 


169.98 


170.30 


169.97 


11.99 


20.54 


127.5975 


0.5837 


87 


0.88 


0.27 


0.022 


3.7244 


-0.322 


0.8208 


4.80 


172.37 


172.43 


172.57 


172.89 


172.56 


12.61 


21.38 


128.3205 


0.5898 


88 


0.88 


0.27 


0.022 


3.7261 


-0.299 


0.8354 


5.80 


167.89 


167.97 


168.09 


168.42 


168.09 


11.59 


19.96 


128.3464 


0.5807 


89 


0.88 


0.27 


0.022 


3.7278 


-0.276 


0.8516 


6.80 


162.92 


163.02 


163.13 


163.47 


163.13 


10.49 


18.42 


128.1998 


0.5695 


90 


0.89 


0.26 


0.022 


3.7244 


-0.297 


0.8439 


6.20 


166.32 


166.41 


166.53 


166.86 


166.53 


11.26 


19.48 


129.1008 


0.5780 


91 


0.89 


0.26 


0.022 


3.7251 


-0.288 


0.8500 


6.60 


164.54 


164.64 


164.76 


165.08 


164.75 


10.86 


18.92 


129.1082 


0.5740 


92 


0.89 


0.26 


0.022 


3.7258 


-0.279 


0.8564 


7.00 


162.70 


162.81 


162.92 


163.25 


162.92 


10.45 


18.35 


129.1188 


0.5695 


93 


0.89 


0.27 


0.022 


3.7257 


-0.320 


0.8171 


3.65 


174.52 


174.55 


174.70 


175.02 


174.69 


13.53 


22.63 


129.0271 


0.5979 


94 


0.89 


0.27 


0.022 


3.7267 


-0.308 


0.8248 


4.20 


172.10 


172.14 


172.28 


172.60 


172.28 


12.97 


21.86 


129.0500 


0.5933 


95 


0.89 


0.27 


0.022 


3.7277 


-0.295 


0.8336 


4.80 


169.36 


169.42 


169.55 


169.87 


169.55 


12.35 


20.99 


129.0430 


0.5884 


96 


0.90 


0.26 


0.022 


3.7245 


-0.312 


0.8294 


4.30 


171.64 


171.68 


171.82 


172.14 


171.82 


12.90 


21.73 


129.7836 


0.5936 


97 


0.90 


0.26 


0.022 


3.7262 


-0.291 


0.8431 


5.25 


167.49 


167.55 


167.68 


168.00 


167.68 


11.94 


20.41 


129.8075 


0.5850 


98 


0.90 


0.26 


0.022 


3.7277 


-0.269 


0.8581 


6.20 


163.12 


163.21 


163.33 


163.65 


163.33 


10.96 


19.04 


129.8293 


0.5756 


99 


0.90 


0.27 


0.022 


3.7274 


-0.314 


0.8164 


2.75 


175.74 


175.76 


175.91 


176.23 


175.91 


14.23 


23.59 


129.7612 


0.6032 


100.... 


0.90 


0.27 


0.022 


3.7276 


-0.312 


0.8177 


2.85 


175.32 


175.34 


175.49 


175.81 


175.49 


14.14 


23.45 


129.7608 


0.6030 


101.... 


0.90 


0.27 


0.022 


3.7277 


-0.311 


0.8183 


2.90 


175.10 


175.12 


175.27 


175.59 


175.27 


14.09 


23.39 


129.7408 


0.6023 


102 


91 


26 


022 


3 725 1 


-0 320 


8194 


2 70 


175 76 


175 78 


175 93 


176 25 


175 93 


14 27 


23 62 


1 30 4990 


6041 


103.... 


0.91 


0.26 


0.022 


3.7264 


-0.304 


0.8302 


3.50 


172.34 


172.36 


172.51 


172.83 


172.51 


13.48 


22.51 


130.4933 


0.5988 


104.... 


0.91 


0.26 


0.022 


3.7277 


-0.286 


0.8418 


4.30 


168.79 


168.83 


168.97 


169.29 


168.97 


12.66 


21.38 


130.5037 


0.5921 


105.... 


0.92 


0.26 


0.022 


3.7271 


-0.309 


0.8222 


2.10 


175.80 


175.81 


175.97 


176.28 


175.97 


14.69 


24.18 


131.1912 


0.6075 


106.... 


0.92 


0.26 


0.022 


3.7274 


-0.306 


0.8243 


2.25 


175.16 


175.17 


175.32 


175.64 


175.32 


14.54 


23.97 


131.2077 


0.6066 


107.... 


0.92 


0.26 


0.022 


3.7277 


-0.303 


0.8263 


2.40 


174.51 


174.52 


174.67 


174.99 


174.67 


14.39 


23.76 


131.1973 


0.6056 


108.... 


0.85 


0.29 


0.023 


3.7258 


-0.320 


0.8169 


6.25 


170.64 


170.72 


170.84 


171.17 


170.84 


11.41 


19.79 


126.1371 


0.5766 


109.... 


0.85 


0.29 


0.023 


3.7265 


-0.312 


0.8223 


6.60 


168.97 


169.06 


169.17 


169.51 


169.18 


11.04 


19.28 


126.1521 


0.5726 


110.... 


0.85 


0.29 


0.023 


3.7271 


-0.302 


0.8288 


7.00 


167.00 


167.11 


167.21 


167.55 


167.22 


10.61 


18.68 


126.1720 


0.5680 


111.... 


0.86 


0.28 


0.023 


3.7244 


-0.314 


0.8276 


6.80 


168.35 


168.44 


168.55 


168.89 


168.56 


10.93 


19.10 


126.9069 


0.5723 


112.... 


0.86 


0.28 


0.023 


3.7246 


-0.312 


0.8291 


6.90 


167.88 


167.98 


168.09 


168.42 


168.09 


10.83 


18.96 


126.8973 


0.5712 
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Table 3 
-Continued 



Model 


M 

(M e ) 


Yi 


z« 


log TV/ 


logi 

(.Lb) 


R 

(Ks) 


Age 
(Gyr) 


(Av > 


(AvO 


<Av 2 > 


<Av 3 > 


<Ay> 


(5v 02 ) 


(■5vi3> 


<Av r > 


r 01 


113.... 


0.86 


0.28 


0.023 


3.7248 


-0.310 


0.8307 


7.00 


167.42 


167.52 


167.63 


167.96 


167.63 


10.72 


18.82 


126.9165 


0.5696 


114.... 


0.87 


0.28 


0.023 


3.7258 


-0.314 


0.8228 


5.55 


170.78 


170.85 


170.98 


171.30 


170.98 


11.91 


20.44 


127.6106 


0.5827 


115.... 


0.87 


0.28 


0.023 


3.7268 


-0.300 


0.8319 


6.15 


168.00 


168.08 


168.20 


168.52 


168.20 


11.29 


19.57 


127.6241 


0.5769 


116.... 


0.87 


0.28 


0.023 


3.7278 


-0.286 


0.8417 


6.75 


165.10 


165.20 


165.31 


165.64 


165.31 


10.65 


18.68 


127.6542 


0.5701 


117.... 


0.88 


0.27 


0.023 


3.7246 


-0.305 


0.8358 


6.25 


167.81 


167.89 


168.01 


168.34 


168.01 


11.27 


19.52 


128.3774 


0.5774 


118.... 


0.88 


0.27 


0.023 


3.7252 


-0.297 


0.8411 


6.60 


166.24 


166.33 


166.44 


166.77 


166.45 


10.91 


19.03 


128.3971 


0.5733 


119.... 


0.88 


0.27 


0.023 


3.7258 


-0.288 


0.8475 


7.00 


164.37 


164.47 


164.58 


164.91 


164.59 


10.50 


18.46 


128.4142 


0.5688 


120.... 


0.89 


0.27 


0.023 


3.7257 


-0.307 


0.8299 


4.95 


170.55 


170.60 


170.74 


171.06 


170.74 


12.32 


20.96 


129.0845 


0.5878 


121.... 


0.89 


0.27 


0.023 


3.7267 


-0.293 


0.8388 


5.55 


167.83 


167.90 


168.03 


168.35 


168.03 


11.70 


20.10 


129.0846 


0.5821 


122.... 


0.89 


0.27 


0.023 


3.7277 


-0.280 


0.8484 


6.15 


165.02 


165.11 


165.22 


165.54 


165.22 


11.07 


19.22 


129.1111 


0.5760 


123.... 


0.90 


0.26 


0.023 


3.7246 


-0.297 


0.8436 


5.70 


167.37 


167.44 


167.57 


167.89 


167.57 


11.62 


19.96 


129.8378 


0.5821 


124.... 


0.90 


0.26 


0.023 


3.7257 


-0.282 


0.8535 


6.35 


164.49 


164.57 


164.69 


165.01 


164.69 


10.96 


19.05 


129.8591 


0.5753 


125.... 


0.90 


0.26 


0.023 


3.7267 


-0 269 


08633 


6.95 


161.71 


16181 


161.92 


162.25 


161.92 


10.34 


18.18 


129.8802 


0.5688 


126 


0.90 


0.26 


0.023 


3.7268 


-0.267 


08641 


7.00 


161.47 


16157 


161.68 


162 01 


161.68 


10.29 


18.11 


129 8680 


0.5682 


127.... 


0.91 


0.26 


0.023 


3.7257 


-0.298 


0.8381 


4.45 


169.93 


169.97 


170.11 


170.42 


170.11 


12.62 


21.33 


130.5190 


0.5917 


128.... 


0.91 


0.26 


0.023 


3.7267 


-0.285 


0.8470 


5.05 


167.28 


167.34 


167.47 


167.78 


167.47 


12.01 


20.49 


130.5456 


0.5861 


129.... 


0.91 


0.26 


0.023 


3.7277 


-0.271 


0.8564 


5.65 


164.54 


164.61 


164.73 


165.05 


164.73 


11.39 


19.62 


130.5533 


0.5805 



Note. — The above mean large separations were calculated averaging over n= 10, 11, 22, 30. The mean small separations were averages over n=10, 11, 12, 30 at 
a fixed / (as indicated). 
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